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Abstract 

To evaluate Riemann's zeta function is important for many investigations related to the area of 
number theory, and to have quickly converging series at hand in particular. We investigate a 
class of summation formulae and find, as a special case, a new proof of a rapidly converging 
series for the Riemann zeta function. The series converges in the entire complex plane, its rate 
of convergence being significantly faster than comparable representations, and so is a useful 
basis for evaluation algorithms. The evaluation of corresponding coefficients is not problematic, 
and precise convergence rates are elaborated in detail. The globally converging series obtained 
allow to reduce Riemann's hypothesis to similar properties on polynomials. And interestingly, 
Laguerre's polynomials form a kind of leitmotif through all sections. 

Keywords: Riemann Zeta function, Kummer function, Laguerre polynomials, Fourier 
transform, Riemann hypothesis. 



1. Introduction and Definitions 

1.1. Confluent Hypergeometric Functions 

Confluent hypergeometric functions are special hypergeometric functions, sometimes called 
also Rummer's function of first and second kind. They are linear independent solutions of Rum- 
mer's differential equation 

zy" (z) + (b-z)y'(z)-ay(z) = 0. 

The first solution, Rummer's function of the first kind, is usually given as a globally converging 
power series 

i=0 [ i ) 

(Rummer's function of the first kind), and Rummer's function of the second kind is often given 
by 

U b >Z> := Jn^'u^ b > $ ~ r( r^ 1) z '~* M (fl-b + l;2-b;z) (1 -2) 

r(l - b + a) T(a) 
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1 INTRODUCTION AND DEFINITIONS 



1.2 Laguerre's Polynomials 



Kummer's transformation states that M(a;b;z) — e z M (b - a;b;—z), which subsequently 
leads to the identity U (a; b; z) - z l ~ h U (1 + a - b; 2 - b;z). 

Another solution of this differential equation - which turns out to be identical to U and thus 
is simply another representation of U - is obtained as an integral representation by the Laplace 
transform 



U(a;b;z) = I ( I +/)""""' d/ 



F(a) 



\-b poo t a-b 

I „-zt 



;dt. (1.3) 



T(l+a-b)J (l + tf 

Some particular and frequently used confluent hypergeometric functions are the upper and 
lower incomplete Gamma function, which have the representations 

= e- z U(l - 5,1 - s,z) = e~ z z s U(l, 1 + s,z) (1.4) 

and 

(-1)* z s+k 



r(*,z):= [ Z f- 1 e- t dt=y 
Jo 



= -M(s,s+l,-z) = -e~ z M(l,s+l,z)- (1-5) 
s s 

1.2. Laguerre's Polynomials 

For a a negative integer the defining series for M reduces to a polynomial, which turns out to 
be closely related to Laguerre's polynomials: The explicit representation is 

LfHz) = ('7)mh> + i,z) -x<-ir(;:*;)|. (i.6) 

In view of (1.2) Laguerre's polynomial may be given by Kummer's second function as well, 
that is L^iz) = tpt/ (-/,«+ l,z). This somehow suggests that Laguerre's polynomial are 
somewhat in between of both solutions M and U of Kummer's differential equation. This is 
central in our investigations and reflected in the results of the next sections. 

In addition to that it is well-know that Laguerre's polynomials are orthogonal with respect to 
the weight-function z a e~ z ; more precisely we find that 

Laguerre's polynomials thus are orthogonal with respect to the inner product 

(g\f)-= / r ? \ u 8(z)f(z)dz. 
Jo r(a+l) 

As a result of the classical theory on Hilbert spaces we may expand a function / in a series 
with respect to this orthogonal basis, provided that ||/|| 2 := \J(f\ f) < oo. The function has the 
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2 FOURIER SERIES 



expansion f(z) — 2i=o/ (ff) ^- a) (z) ( an d is therefore in the closed hull of all L w ), its coefficients 
taking the explicit form 

r L (a) (z) z a e - z 

(cf. (1.7)). Conversely, the norm can be recovered from the function's coefficients, as \\f\\\ = 
An elementary example of such a representation in explicit terms is 



(1 



as can be verified straight forward by evaluating the respective integrals for the coefficients (1.8). 
Notably, this series converges point-wise if t > — | (H/H2 < °°), even if a < -1. 

To give a reference of these classical ingredients aggregated in this section above we would 
like to refer to the standard work [AS64]. 



2. Fourier Series 



Kummer's functions allow some explicit representation as series of Laguerre polynomials. 



Theorem 1 (Expansion of Kummer's functions in terms of Laguerre polynomials). Suppose that 

5 

~v 



Re (b - a) > and Re (a - 2b) > — 5, then 



tb-l\ t-a\ 



and 



moreover, 



/or Re (.-§)< 



nt u \ n (,),. (i) . 
U(a;b;z) = - - rr / , ^; fe) 77 ~ (2-1) 

(1 + a - b + a)\ j-£ ' fb-a-a-2\ 



r(,; Z ) = zV*Y ^' U ; - (2.2) 



Proof. As for the proof notice first that ^ = (-1)' 2}=o (j-"}^f~^ wn i cn I s a °f con " 
verse of (1.6) and verified straight forward by comparing the coefficients of respective powers 
of z. Substituting this into (1.1), interchanging the order of summation and employing Gauss' 

3 



2 FOURIER SERIES 



hypergeometric theorem (cf. [Koe98]) gives 



M(a;b;z) = ^ 



(" 



■(-.yz(,?>r« 

oo oo /fl+i-l\ . , 



(=0 i 



2 l<M ( Z ) r w r (fr -«-/?)!> + ;) 

7=0 

J>^( Z )(-1) 



./=<> 



(V) 



(b-p-\\ 

\ a ) J=0 



T{d)T{b-a)T{b-p + j) 

, n m 
mm 

(7) 



which is the desired assertion. As regards convergence it follows from the elementary recursion 



(2.3) 



that (- 1 V L J ;) (z) - Oi^j 13 1 ), a more thorough analysis even discloses that - = e z (l +0(j 

Hence, Lj ^ (z) — o{j~P~ l j^Fi") = ^( ]F5t)> anc ' tne ser i es converges - irrespective 
ofyS-ifRe(fc-a) >0. ' 

To verify the second statement notice that 



U? = 



U (a; b; z) dz 



Lf\z) z a e-- 
T(a+ 1) 

Lf\z) z a e-i 1 

o r(a + i)r( fl )y 



-- z >f- l (l + t) b - a - l dtdz 



o r («) 
i 



r(a) 7 



f"- 1 (1 + f) 



(1 + 0' 



—df, 



Q-+I+ 1 



which is a consequence of (1.3) and (1.9). The latter integral is a beta function, we thus continue 
and find the coefficient 
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2 FOURIER SERIES 



Algorithm 1 to evaluate Lf ] (z), based on Lf (z) = (2 + (z) - (l + ^f 1 )^ £>• 

L1:=0; Laguerre:= 1 
For j : = 1 to i 

L0:= LI; Ll:= Laguerre ; 

Laguerre: = ((2* j+ alpha- 1- z)* LI- (j+ alpha- 1)* LO)/ j; 

Next j 

Return Laguerre 



Algorithm 2 to evaluate hf l) (z), based on (2.3) 

L1:=0; Laguerre:= 1 
For j : = 1 to i 

L0:= LI; Ll:= Laguerre; 

Laguerre: = ((beta+ 1- j- z)* LI- z* LO)/ j; 

Next j 

Return Laguerre 



Uf = -J— I f~ 1+i (1 + tf-"-"-'- 2 df 



Y(a) 

1 Y(a + i)Y{a + 2-b) 
r (a) F (i + a + a + 2 - b) 

(a+i-l\ 

= {a+\-b)\ { j ) 
(i + a + a + 1 - b)l ft+a+<*+2-b\ 

which is the respective coefficient for U. 

Convergence is more difficult compared to Rummer's function of the first kind. However, 

L (a \z) 

we will show below (Theorem 9) that lim sup, * g t < oo and the series thus convergences if 
Re(a-2fc) > -5. □ 



Remark 2. Identity (2.1) can be found for the special case a = b — 1 in [EMOT53]. 

To evaluate the series above using Laguerre polynomials it is necessary to have good eval- 
uations of Laguerre polynomials at hand for given parameters a and z. Moreover, numerical 
approximations should be sufficiently good and the computation stable. Although there are ex- 
plicit expressions or Horner's scheme available to evaluate Laguerre polynomials, the resulting 
algorithms usually behave unstable very soon. 

We have found the algorithms described much better, adaptations even allow to evaluate and 
then successively store intermediary results. However, they should not be interchanged, as this 
will cause numerical instability again. 
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3 CONTINUOUS FOURIER TRANSFORM AND POISSON SUMMATION FORMULA 



3. Continuous Fourier Transform and Poisson Summation Formula 

3.1. Continuous Fourier Transform. 

It is well-known that Poisson's summation formula provides an efficient tool to evaluate 
sums, some authors dedicate entire chapters to these summation techniques, see for instance 
[IK04]. To apply these effective summation identities we need to have a good expression for the 
functions involved at hand, which involve the continuous Fourier transform. 

In literature there occur a few variants for the continuous Fourier transform of a function /, 
which are interchanged frequently; for our purposes it is most convenient to state 



/oo 
f(z)e- 2nik Az 
oo 



as a definition. 

We will sometimes abuse the notation just introduced and improperly write T f (z) shortly 
for T (f (z)) (z), that is to say we use the argument z for both functions - / and T (/) - synony- 
mously. 

The following result establishes that Laguerre's polynomials, as well as Kummer's functions 
are each others Fourier transform in the following sense: 

Theorem 3 (Fourier transform of Kummer's functions). 

1. Te-^Lf (z 2 n) = {-Ije-^L)-™-^ (z 2 n), 

2. Te- z2 "u(a,b,z 2 n) = ^k^e^M (a,a + § - b,z 2 n), 

3. T e- £ -"M[a,b,z 2 n) = ^^e~ en u[a,a+ § - b,z 2 n). 

Proof. Notice first, that e~ z ~ K is an eigenfunction of the Fourier transform T, which already cov- 
ers the desired statement for i = 0. Now recall that - due to integration by parts - T ( £jf (zfj = 
(2niz) J T (/ (z)) and observe that differentiating this first eigenfunction involves Laguerre's poly- 

>2j 2 _ 2 (— — ) / 9 \ 

nomials again, as z " = {-An) 1 j\e z n L\ 2 \z nj. Hence, 

(z 2 n) J 2 (_!); 2 t_iy I d 2j \ 

y -^-e- rn = -M- {2niz) 2j T = e ' — — '/ \-^r : r \ 
]\ {An) J j\ (An) 1 j\ \dz 2 J / 

- S7j! r h-yfl^ (A)) - r<*$ " (*)■ 

Next, by (1.6) and the identity just derived, 



(-1)' , ' (z 2 n) = <** g {-if ^ 3 j ! 



=r^^-j:;-^(fn) = r^Lf^n), 
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4 APPLICATION TO RIEMANN ZETA FUNCTION 



3.2 Poisson Summation Formula. 



the latter identity being a special case (jc = 0,/3= -|) of the more general identity 

J>g(;t)Lf (y) = &^\ X + y). 

This proves the first statement. 

The other statements are an immediate consequence of this first one and Theorem 1 . □ 

3.2. Poisson Summation Formula. 

Given the notation introduced Poisson's summation formula reads Zzez/fe) = Tikezf(k) 
(cf. [Zyg68, Pin02]), or a bit more generally 

Y J f((z-zo)x)e-^-^ = l£/^±M e -**o (3>1) 

zeZ X kel \ X ' 

when applied to the same function / but translated (zo), contracted (x) and shifted in phase (ko). 

As already mentioned, the Poisson summation formula often transforms slowly converging 
series into very rapidly converging series. We will exploit this fact to obtain the following results, 
which are an application of Poisson's summation formula to the Fourier transforms elaborated in 
the latter section. 



4. Application to Riemann Zeta Function 

It turns out that a particular application of summation techniques outlined above is a very 
rapidly converging series for the Riemann zeta function. We are even able to prove these follow- 
ing variants, and Riemann's function equation follows as a by-product: 

Theorem 4. For any s e C we have 

s(s-l) ^ = 1-5(1-^ V - / + V ,., ' . (4.1) 

More generally, for every arbitrarily chosen x satisfying Re (x) > |Im (x)\ we find 



= V 1 - S) X + S X 

r( f ,tt) r(!=*,fF) 



s(s-l) '- =(1 - s)x s + sx s ~ l + 



k=l 



moreover 



(2 s - 1) (l - 2 S J = T x (s) + Ti (1 - s) , 



7H 



where T v (s) is the 2 forward difference 

Y^s {Ak+\fn£ \ p^s (4k+2) 2 nx 2 j p^s (4k+3) 2 nx 2 j 



A=( , ^ (4t+l) 2 ;r j 2 ^ (4fr+2) 2 ;r j ; ^ (4fr+3) 2 ;r j 2 
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4 APPLICATION TO RIEMANN ZETA FUNCTION 



Remark 5. It should be stressed that these are globally convergent series for a function analytic in 
the entire plane, converging in particular for s = 1. Moreover, by de l'Hopital's rule, lim -tt-% = 

lim ._, >e ,~. - = 1, the rate of convergence thus is of order 

7 s '? z— (s—l)z s -e - ' & 



(CO \ 

JV^ = 0(e- k2 ») 
k'=k ' 

= (9(0.04321 ... k2 ), 



which is pretty quick. 

As an additional result, Riemann's function equation follows immediately from equation 
(4.1). 

It seems that an identity close to (4.1) already was known to Riemann himself, although the 
proof being based on integrals involving the Jacobi function #3 (z) := TikeZ e~ k " z , or rather the 
function 6(z) := (z) + \ ff{ (z): Both satisfy 6 3 (z) = ^6» 3 (~), but (k) (z) -> as z -> and 
^ — > 00 for any A: e {0, 1, 2, . . . }. 1 By exploiting this fact (4.1) can be recovered as well - this 
is for example demonstrated in the third method (out of 7) of proving the function equation in 
[Tit86]. 

Remark 6. The famous, fast algorithm for multiple evaluations of the Riemann Zeta function 
in [OS88] is based on an a representation for the zeta function, which is not converging on the 
entire complex plane, and the rate of convergence being slower, but allowing precise bounds and 
estimates. t 

The representation given by [Son94], £(s) = H„=i ZLo Wifw' conver 8 es globally 

as well, however, the rate of convergence is significantly slower. 

A similar algorithm is proposed in [Bor91] with exponential convergence rate. 

So to summarize and compare the representations described in Theorem 4 are faster, the only 
downside is that r(f) is very small for huge imaginary parts, and incomplete Gamma functions 
have to be evaluated. However, this can be done very quickly, as will be further outlined below. 

Proof. As for the proof notice first that 

f e -^M(a,b,^)= 1 --^- f e-^u(a,a + 3 --b,^], 
t-i v 1 x T (b - a) t-L \ 2 x 2 I 

Z— — DO J— — CO ' ' 

or, using (1.2), 

1 + 2 1 e~ M M (a, b, z 2 ttx 2 ) 

-\m^i^^7^u{a,a + l-b^)\ 

We choose a — 1 and b = 1 + I. Thus, using (1.5) and (1.4), 

1 + 2 s_f r(|A 2 ) = s + 2 r(i + f) ^ r(^, % 



(z 2 7Tx 2 y (s-i)x x r (|) j-* 



'By the way: Suppose that 9{z) = Z'"' 1 • then (-1)' ^p-S(z) = £} =0 (ty)^^ (i)> strikingly reminding 
to the explicit representation of Laguerre polynomials (1.6). 



4 APPLICATION TO RIEMANN ZETA FUNCTION 



or 



- . - ,y 4H y r (¥*») 



This nice identity at hand we find that 



((s)TU) yU-k 2 nx 2 ) + TU-k 2 nx 2 ) 

s(s-i) — ^i=i(j-i)y 



(1 -s)x s + sx s ~ 1 



si 1 " 1 



which is the second statement. Notice, that the condition Re (x) > |Im (x)| insures both, Re (x 2 ) > 
and Re \X\ > 0, which is necessary for convergence. 

The first statement is obvious by the choice x — 1 . 

As for the next statement recall that 

1 ' - H (k 2 nY 

In order to get rid of slowly converging y (and replace it by rapidly converging Fourier transform 
T) we apply Poisson summation formula (3.1) again, now with zo = \ and ko = 0. Hence, 

{ 1 * f 5 -(M) 2 -) ¥ ( ' 27r)f 

il 1 y r (^^) |(1 /(f^) 

* 2 ' 5 2 ' S "^((2 fe -l) 2 ,)- W ! ' 

as the statement holds for | as well. Combining the latter two identities and rearranging the 
terms is cumbersome, but finally gives the assertion. □ 



Remark 7. x := -§y is notably a possible choice if Im(s) > Mp 2 - ~ 1.21, as in this case 
Re (x) > |Im (jc)|. The advantage of this this particular choice is that the term (1 - s) x s + s x^ 1 
vanishes. 

Remark 8. It should be noticed that the identity (4.2) is central here. It allows to replace the 
slowly converging series (which involves y) by its Fourier transform, which is the very fast 
converging series (which involves F). This is the key strategy in all improvements of convergence 
above. 
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5 ASYMPTOTICS OF THE LAGUERRE POLYNOMIALS 



5. Asymptotics of the Laguerre Polynomials 

Theorem 9. (Asymptotics of the Laguerre Polynomial) Let Re (z) > 0. Then, as i — > oo, 
L f 6) * ^ j!r cos (2 ^(Z+V) " ! (« + I)) «** 

4 a) (-^#^*p(2#+^))- 

Remark 10. We write /j « g,- as / — > 00 to indicate that lim/^oo ^ = 1. 

Remark 11. To prove the statements of the theorem we will involve Bessel functions. Bessel 
functions, of first and second kind, are 

(-1)"' lx\ a+2m 1 (x\ a+2m 

J a (x) := V ; and7 a x) : =V— — . 

^ mlT(a + m + 1) \2I m lT(a + m + 1) Ml 

Notice, that = are entire, even functions. Moreover, Kummer's second formula links 

\IX) XT 

Bessel functions to confluent hypergeometric functions, as j^y = e~ : ^"^(a+x) 2 ^ • 

Proof. We start with identity 13.3.7 from [AS64] which states that 

M{a,b,z) z lz\" Jb-Un(i2z(b-2a)) 

— — — = e> ■ 2^A n ■ I- I , (5.D 

1 {b) 7~t) yz/ Q V2z(fc-2a)) 

where 7 is the Bessel function of the first kind and A satisfies the recursion Aq := l,Aj := 0, A2 := 
I and A„ := (l + ^) A„_ 2 - ^A„_ 3 . Although we refer to this Identity (5.1) without proof 
we mention that it follows straight forward by successively comparing the respective coefficients 
of z in (5.1). The coefficients satisfy A„ = 0\r~ b j\ for b > 2a, as follows directly from the 
recursive definition. 
Following [AW05], 

J a (x) « , /— cos (x - jj- (a + I- 1) , I a (x) x (5.2) 



2// y/2nx 
as x — > 00. 

Combining all those ingredients we see that the initial term (n = 0) dominates all other 
summands in (5.1) when a — > -00, so we may neglect them for n > 1 to identify the rate of 
convergence. But for n = the assertion states that 



M{a,b,z) i Jb-i(y/2z(b-2a)) 
r(b) {^2z(b-2a)) 



b-l ' 



which is a useful approximation itself. 

There is another, similar approach, which is quite useful, which we want to give here as well 
for the sake of completeness and further reference: It starts with identity 13.3.8 from [AS64] 
with parameter h — | . This reads 

M(a,b,z) _ z YV n-Jb-i+nfeypaz) 

r(b) - e2 'h CnZ ( . — ' 
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6 APPLICATIONS TO THE INCOMPLETE GAMMA FUNCTION 



where C = 1, d = -b, C 2 = ^ and C„ = -^C„_! + (I + ^)C„_ 2 + ^C„_ 3 . Here, 
C„ = O ^{nfj and similar to the reasoning above we obtain 

M(a,b,z) : Jb-i(2^Paz) 

From both asymptotic identities the theorem follows in view of (1.6), that is 

L< B) fe):=('' + a )jlf(Hff+U), 

and (5.2) as z — > oo. □ 

Remark 12. We want to stress that the method used in the proof above allows to compute terms 
of higher order of the expressions given in Theorem 9. The higher order correction terms give 
successive improvements of order 4=. 

For an interesting treatment to evaluate Laguerre polynomials for large i we refer the reader 
to[BBC08]. 



6. Applications to the incomplete Gamma Function 

In order to make use of the rapidly converging series (4.1) it is necessary to have a good 
implementation for the upper incomplete gamma function at hand. To this end we further elabo- 
rate on a procedure based on continued fractions (cf. [AT99]), which is always a good candidate 
for rapid convergence. The formula is a special case of Gauss' continued fraction method using 
confluent hypergeometric functions (see [JT80, Wal48]): 

z s e- z 

r(*,z) = (6.1) 

1 — s 

z + 



2-s 



1 + 



2 



3-s 

z + 

3 

1 + 



z+ ■■ 

(as aformal power series this is equivalent to T(s,z) = z s e~ z ■ 2-F0 (1, 1-5, - jj). 
Stopping a general continued fraction gives its convergents, the k' h -convergent is 

Pk fl i 



Ik a 2 
bi + 



b 2 + 



Ok 

11 h 



6 APPLICATIONS TO THE INCOMPLETE GAMMA FUNCTION 



To simplify the respective values there is - from standard theory [Wal48] - the recursion 

ft W Pk - 2 ) + h( ^ 



qk I \ qk-i I \ qk-\ 



with initial conditions ( ^ I := ( ? I and [ Pl I = ( f 1 |, and additionally the somewhat 
\q ) [1 ) \qt ) \b! } 

more explicit formula 



p* -y (-iy-1 " (6 . 2) 



?* £f q> ■ qi-\ 



When applied to the continuous fraction of the incomplete gamma function (6.1) the recursions 
for the numerator and denominator Pk \ s -~\ read 



Pk 0, z) 

and 



[(§ - sjpk-i{s,z) +Pk-i (s,z) fceven 
[ h Y L Pk-2(s,Z)+Zpk-i(s,z) kodd 

\(%-s)qk-2(s,z) + qk-i{s,z) A: even 



q k (s,z) = \\U * % , ' (6.3) 

but different initial values: po — 0, p\ = 1, q-\ = and go = 1. It comes without surprise that 
Laguerre polynomials appear again, as these recursions have the closed expression 

qik (s, z) = k\L { k s) {-z) and q 2 k+\ (s, z) = k\z Lf~ s) (-z) . 

In view of (6.2) the knowledge of the denominator is already sufficient to get the convergents, 
in explicit terms we get the simple expressions 

P2k(s,Z) ri (7) T(s,z) 

q*M A L (-,) ( _ z)L (-,) ( _ z) z s e -z 

and 

P2M (s,z) = 1 + 1 y (V) > r(^,z) 

to + i(*,z) z ztiI^HiP _,) (-?) 

which have been observed in [BBC08]. To find a handy expression for the numerator as well is 
surprisingly much more difficult - but yes, the somewhat curious result involves Laguerres again: 

Theorem 13. For any s and z + we have 

y L¥J /( 2i- S ) , , (Tl) 

T(s,z) ,. Li = l ^k-li+l^ z > 
= l im 

V L¥J jVm-s) r _^C£) 

1 A'=l H-2i+\ ^ Z > (t\ 

= - + Hm jj-, (6.5) 

g yg 1. 

12 



7 APPLICATION TO RIEMANN'S ZETA FUNCTION 



Remark. The integer valued floor function [.J satisfies x — 1 < L*J ^ x. 

Proof. The fractions in the limit represent explicit expressions for ( gzTifo'zj ' res P-)- Again, 

the proof links the recursions (6.3) introduced above to well known recursions of Laguerre poly- 
nomials. 

The rate of approximation is an interesting consequence of (6.2) (same for (6.4)): 



O 



Z\j T (-s. 



1 

i+1 



(7) 



\i=k 



= o 



o 



' 00 — J-l 

5 rHiW^ 



d/ 



o(r(o,4Vfe)) 



= 



which is finally the desired rate. 



□ 



7. Application to Riemann's Zeta Function 

We have given a few approximations for the upper incomplete Gamma function which in- 
volve polynomials, for example (2.2), (6.2) ((6.4), respectively) and (6.5). Those expressions 
converge sufficiently quick and do not impose any difficulties as their argument z tends to infin- 
ity. So they can be substituted in (4.1) to give recent approximations and a variety of possible 
investigations on the geometry of the Riemann zeta function. 

As an initial example combine (2.2) and (4.1) to get the representation 



£ A) (*):= 1 -*(!-*)£ 



n-1 

„(A) 



k=0 



1 1 



\\ k+l ) { jfc+1 )) 



< n f( ' )r (*) 

-> s ( s - 1) - , 



where :- - £ ^— — ^ — ¥Jll t Notice, that this representation converges, as n — > 00, uniformly 
on -2A - 1 < Re (s) < 2 A. 

Moreover - and this is a key observation - the numerator or the function 

^§ + «-l + Ap + «-l + A 

as a function in variable s, is a polynomial of degree 2«. Thus, by Hurwitz' Theorem in complex 
analysis (cf. [Con78]), the zeros of £ are just the accumulation points of s zeros. 

Interestingly, the zeros of the polynomials above have a very nice, symmetric pattern in 
common, a typical result is plotted in figure 7.1 2 : The total of zeros is 50, 18 on the left (right) 



^Figures have been computed using Mathematica. 
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Figure 7.1: Zeros of ft?. 



of Re (s) = 1 7 have a positive (negative) imaginary part on the critical line, a few exemptions 
occur close to 2 A (1 - 2 A, respectively). 

However: all zeros in the area of convergence, which is -2A - 1 < Re (s) < 2 A, lie on the 
critical line. 

Another pattern is found when employing (6.2) or (6.5), in our next example for k — 6 and 
involving z 2 n in a, for z from 1 to 5: The corresponding polynomial has degree 48, its zeros are 
depicted in figure 7.2 (6 do not fit to the scale chosen). 

Interestingly, the most convenient pattern is observed when substituting (2.2) into the equa- 

L A 

tion Yjc (s). Stopping the infinite sum including — t * 1+g _i as a final term and solving Tj (s) + 

Yi (1 — s) = again leads to finding the roots of polynomials in s, which a of degree 2k here. Ob- 
viously, the zeros corresponding to the factor (2 s - 1) (l - 2'~ s ) have to become visible as well, 
but, as numerical experiments show, those zeros are being added very slowly, as k increases. 
We have depicted some zeros for k — 50 in figure 7.3, where only 4 zeros corresponding to 
s — 2±3±^§(^3 ~ 9.06) appear. However, and this is potentially a big advantage in com- 
parison to the other figures, all other figures are located on the critical line and two symmetric 
bubbles except a few others usually located close to the real line. The complete pattern of zeros 
is depicted in figure 7.4. Do these associated polynomials always have their zeros in the region 
of convergence on the critical strip? The answer "Yes" obviously implies RH. So this ques- 
tion seems being worth an attempt and has to be investigated in much more detail in promising, 
further research. 
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Figure 7.4: Zeros of a = 5 
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